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Abstract 


The  equilibrium  geometries  and  ionization  potentials  have 
been  calculated  self-consistently  for  lithium  clusters  con¬ 
sisting  of  up  to  five  atoms.  The  exchange  interaction  is  treated 
exactly  while  the  correlation  effects  are  included  perturbatively- 
The  correlation  contribution  is  found  to  have  important  influence 
on  the  magnetism  of  the  cluster  and  its  relationship  with  the 
topology. 


There  has  been  a  growing  interest  in  the  study  of  the 
electronic  properties  of  small  metal  clusters  in  recent  years. 

This  is  prompted  not  only  because  small  metal  particles  may 
have  technological  importance  but  also  because  their  study  is 
interdisciplinary  and  bridges  the  gap  in  our  knowledge  of  the 
electronic  structures  between  atoms  and  solid  state  systems. 

Recent  improvement  in  the  experimental  techniques  ■*  has  enabled 
one  to  produce  clusters  of  atoms  ranging  from  2  to  100  atoms  per 
cluster.  Three  important  features  are  noticed.  (1)  Conspicuously 
large  peaks  occur  in  the  mass  spectra  of  the  clusters  indicating 
that  certain  clusters  of  atoms  are  more  abundant  than  others. 

(2)  These  "magic  numbers"  (i.e.  the  number  of  atoms  in  the  more 
abundant  clusters)  are  dependent  not  only  on  the  type  of  atoms 
but  also  on  experimental  conditions.  (3)  The  clusters  con¬ 
taining  even  number  of  atoms  are  more  abundant  than  those  con¬ 
taining  odd-number  of  atoms. 

In  a  recent  experiment  on  sodium  clusters,  Knight  et.  al.^ 
found  the  magic  numbers  for  N  =  2,  8,  20,  40,  58,  and  92  where 
N  is  the  number  of  atoms  in  the  cluster.  While  these  results 

are  somewhat  different  from  that  published  earlier  by  Kappes 
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et.  al.  ,  the  remarkable  aspect  of  the  work  of  Knight  et.  al. 

is  that  they  could  explain  the  occurrence  of  these  magic  numbers 

in  terms  of  a  simple  one-electron  shell  model  based  on  jellium 

calculations.  This  implies  that  the  occurrence  of  above  magic 

numbers  should  be  a  universal  feature  of  all  simple  monovalent 


metals. 


In  this  letter,  we  report  three  different  calculations  for 
lithium  clusters.  (1)  We  have  carried  out  spin-polarized  jellium 
cluster  calculations  based  on  local  density  approximation  for 
N  =  1  to  42.  The  second  derivative  of  the  total  energy  showed 
peaks  at  N  =  2,8,20,  and  40  in  agreement  with  the  calculations  of 
Knight  et.  al.^  These  correspond  to  clusters  with  filled  electron 
shells.  In  addition,  we  also  observed  additional  peaks  at  N  =  5, 

13,  27,  and  37  that  are  absent  in  the  calculation  of  Knight,  et.  al. 
These  peaks  originate  from  the  half-filled  shells  and  the  filling  of 
the  orbital  spin  states  follow  Hund’s  rule.  Thus,  if  the  jellium 
model  is  correct,  the  results  would  imply  that  N  =  13  cluster  should 
be  more  abundant  than  the  12  and  the  14  atom  clusters.  The  experi¬ 
mental  result^  is  just  the  opposite.  (2)  To  resolve  this  issue  and 
to  try  to  understand  the  success  of  the  original  jellium  model, 
we  have  carried  out  unrestricted  Hartree-Fock  (UHF)  calculations 
for  Li  clusters  consisting  of  one  to  five  atoms.  The  total  energy/ 
atom  decreases  continuously  as  N  increases  from  1  to  5,  thus 
revealing  no  tendency  for  N  s  2  to  be  a  magic  number.  When  inter¬ 
atomic  correlation  was  included,  the  total  energy/atom  for  N  =  2 
did  exhibit  a  dip  compared  to  N  =  1  or  N  =  3.  This  implies  that 
correlation  effects  are  crucial  in  studying  the  relative  binding 
energy  of  small  metal  clusters.  The  equilibrium  geometries  of 
all  our  clusters  in  the  calculations  containing  correlation  are 
planar,  in  agreement  with  other  observations^  These  low 
symmetry  planar  geometries  are  a  result  of  the  Jahn-Teller  effect 
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which  connects,  via  the  Hellmann-Feynman  theorem^,  the  equili¬ 
brium  nuclear  positions  with  the  electron  density  provided  by 
the  molecular  orbitals.  For  the  4-atom  clusters,  we  have  con¬ 
sidered  two  possible  spin  configurations,  singlet  and  triplet. 

Vie  find  that  the  ground  state  is  nonmagnetic  (spin  singlet)  . 

This  is  in  contrast  with  the  jellium  calculation  which  predicts 
the  triplet  state  to  have  lower  energy.  The  source  of  this  dis¬ 
crepancy  is  found  to  lie  in  the  nature  of  preferred  filling  of 
the  molecular  orbitals  as  the  cluster  size  increases. (3)  The 
ionization  potentials  of  these  clusters  have  also  been  calcu¬ 
lated.  The  UKF  results  for  these  decrease  monotonically  with 
the  increase  of  cluster  size.  When  correlation  is  included, 

g 

structures  similar  to  that  observed  in  sodium  clusters  appear. 
This  reinforces  our  earlier  claim  that  correlation  effects  are 
essential  in  understanding  the  electronic  properties. 

We  now  present  the  details  of  our  calculation.  The  jellium 
calculations  were  carried  out  by  representing  the  external  per¬ 
turbation  by  a  spherical  distribution  of  homogeneous  density  of 
positive  charge,  namely, 

n  ^  ir)  =  0  (R-?)  ,  (D 

ext  o 

where  R  =  r  is  the  radius  of  the  cluster.  The  electron 

s 

radius  r  and  electron  density  n  are  related  through  the 
s  ^ 

equation  1/n  =  (4Ti/,)r  The  electronic  structure  of  the  cluster 

O  V  s 

of  K-atoms  was  calculated  self-consistently  by  solving  the  den¬ 
sity  functional  equation, 


iliUi.iA.Aii. 


The  Dotential  is  composed  of  electrostatic  and  exchange- 

correlation  potential  for  spin  o  in  the  local  spin  density  approxi¬ 
mation.  and  ij;?  are  respectively  the  energy  eigenvalue  and 
wave  function  of  the  i^^  electron  with  spin  o.  The  total 
energies  of  the  clusters  E(N)  for  1^N$42  were  calculated.  In 

II 

Fig.  1  we  plot  the  second  derivative  of  this  quantity  E  (N)  = 

E (N+1) -2E (N) +E(N-1)  in  order  to  fascilitate  the  comparison  of 
our  spin  polarized  calculation  with  the  result  of  Knight  et.  al.^ 
The  peaks  at  N  =  2,8,20,34,  and  40  correspond  to  the  filling  of 
ls,lslp;  lslpld2s;  lslpld2slf  and  Islpld2slf2p  shells  as  in  the 
work  of  Knight  et.al.^  We  observe  additional  peaks  at  N  =  5, 

13,27,  and  37  which  correspond  to  half-filling  of  orbitals  Ip, 

Id,  If,  and  2p  orbitals  in  keeping  with  Hund's  rule  coupling. 

The  UHF  calculations  were  carried  out  by  representing  the 
molecular  orbitals  as  a  linear  combination  of  the  Gaussian-type 
atomic  orbitals  centered  at  the  atomic  sites.  Interatomic  corre- 
lation  correction  was  introduced  perturbatively  by  including 
pair  excitations  of  the  valence  electrons.  The  equilibrium  geo¬ 
metries  for  both  neutral  and  singly  ionized  clusters  were  obtained 
by  calculating  the  Hellmann-Feynmann  forces  on  every  atom  and 
following  the  path  of  steepest  descent  to  the  minimum  energy  con¬ 
figuration.  The  equilibrium  geometries  are  given  in  Fig.  2,  and  the 
corresponding  bond  lengths  are  listed  in  Table  1.  The  energy/atom 
of  the  neutral  cluster  obtained  from  UHF  as  well  as  with  correlation 
correction  are  given  in  Table  1.  Note  that  in  UHF  calculation 
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the  energy/aton  decreases  continuously  unlike  that  in  the  jellium 
calculation.  However,  when  correlation  is  introduced,  the  energy/ 
atcni  for  N  =  2  cluster  is  distinctly  lower  than  those  for  N  =  1 
and  N  =  3.  Thus,  purely  on  energy  considerations,  one  would  ex- 
nect  a  dimer  to  be  more  abundant  than  either  a  monomer  or  a  trimer. 

This  establishes  N  =  2  as  a  magic  number. 

Using  the  optimized  geometries  of  the  neutral  and  the  ionized 
clusters  the  ionization  potentials  were  calculated.  These  are  pre¬ 
sented  in  Table  1.  The  experimental  ionization  potential^®  of 
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lithium  is  5.39  eV.  For  N  =  2,  the  impact  and  photo-ionization 

aocearance  potentials  are  4.8  eV  and  5.1  eV  respectively.  For  N  =  3, 

the  imoact  appearance  potential  is  4.3  eV.  Experimental  results  for 

higher  clusters  are  not  yet  available.  However,  the  nature  of 

variation  of  the  ionization  potential  with  cluster  size  obtained 

8 

with  correlation  is  consistent  with  that  found  in  Sodium  .  The  equili¬ 
brium  geometries  and  ionization  potentials  for  3  are  in  good  agree¬ 
ment  with  previous  calculations^^.  The  fact  that  the  small  even 
clusters  have  larger  ionization  potentials  than  the  odd  ones  suggests 
that  they  are  less  reactive  than  the  odd  clusters.  This  leads  to  the 

observed  results  that  the  even  clusters  are  more  abundant  than  the 

, ,  1,4 

oad  ones 

Ive  now  discuss  the  magnetism  of  the  clusters.  The  minimum 
value  of  N  where  we  can  make  a  choice  of  spin  populations  is  4.  In 
this  cluster,  one  can  either  have  +  t  (spin  triplet)  or  t  +  +  l 

(spin  singlet)  configuration  for  the  valence  electron  spins. 

Taking  the  planar  cluster  obtained  for  the  spin  singlet  case 
(shown  in  Fig.  2  (iii))  we  reoptimized  the  bond  lengths  and  angles 
for  the  triplet  case.  This  energy  was  0.0164  eV  higher 
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than  the  singlet  energy  (including  correlation)  .  In  the  jelliuin 
calculation,  however,  the  triplet  was  found  to  be  the  ground  state. 

To  see  if  this  effect  is  associated  with  the  dimensionality  of  the 
clusters  (jellium  is  three-dimensional  whereas  Fig.  2  (iii)  is 
planar)  we  calculated  the  energy  for  both  spin  configurations  by 
optimizing  the  bond  lengths  of  the  N=4  cluster  where  the  atoms  are 
arranged  in  a  bcc  tetrahedron.  The  energy/atom  of  the  tetrahedral 
cluster  for  the  spin  triplet  is  -202.4884  eV  while  for  the  spin 
singlet  tetrahedral  cluster,  the  energy/atom  is  -202.0729  eV.  Thus, 
the  three-dimensional  N=4  cluster  prefers  to  be  magnetic.  This  is 
in  agreement  with  the  spherical  jellium  model  calculations.  V?e  must 
remind  the  reader  that  the  planar  configuration  (Fig.  2 (iii)  is  the 
equilibrium  geometry  obtained  by  complete  optimization.  The  pre¬ 
ferred  spin  orientation  and  its  relationship  with  the  dimensionality 
of  the  cluster  can  be  understood  in  the  following  way.  An  analysis 
of  the  populations  of  the  molecular  orbitals  of  the  planar  clusters 
reveals  that  the  p  -,  p  -  and  p  -like  states  are  filled  in  succession. 
For  the  N=4  cluster,  p^-  and  p^-like  states  are  empty.  For  the  three- 

dimensional  tetrahedral  cluster,  on  the  other  hand,  p-,  p-,  and  p-  like 

X  j  z 

states  are  filled  simultaneously  and  hence  energy  can  be  lowered  by 
maximizing  spin  due  to  exchange  interaction.  In  addition,  we  want 
to  point  out  that  the  spin-singlet  N=4  planar  cluster  is  more  com¬ 
pact  than  the  spin-triplet  N=4  tetrahedral  cluster.  This  leads  to 
a  larger  electron  density  along  the  bond  in  the  two-dimensional  clus¬ 
ter  as  compared  to  the  corresponding  three-dimensional  one.  The  in¬ 
creased  electron  density  favors  anti-symmetry  in  the  spin,  making 
the  planar  cluster  less  magnetic.  Since  the  dimensionality  of  the 
cluster  topology  is  linked  to  the  magnetic  order  for  N  <  7  in  mono- 


valent  metals,  we  can  say  that  the  jellium  model  would  predict  the 
wrong  spin  configuration  for  N  -  4,  5,  6  clusters.  We  have  found 
that  the  equilibrium  geometry  for  N  =  5  in  Fig.  2(iv)  has  th'^  spin 
orientation  f-l-  +  instead  of  ++  t  +  t.  As  a  further  test,  we  have 
also  calculated  the  ground  state  geometry  of  the  four  atom  beryl¬ 
lium  cluster.  In  this  case,  there  are  8  valence  electrons  and 
the  complete  filling  of  the  p-like  states  should  result  in  three- 
dimensional  clusters.  As  expected,  we  did  obtain  the  equilibrium 
cluster  as  a  tetrahedral  arrangement  of  the  four  Be  atoms  in¬ 
stead  of  a  planar  one.  The  interplay  between  the  electronic  and 
the  ionic  arrangements  may  also  reduce  the  total  spin  predicted 
by  the  jellium  model  in  larger  clusters  than  the  ones  considered 
here.  This  would  explain  why  the  half  filled  shells  do  not  show 
up  as  magic  numbers.  Another  explanation  might  be  the  reduction 
of  the  formation  probability  of  odd  numbered  clusters  as  a  result 
of  larger  reaction  rates^^  of  odd  clusters  as  compared  to  the  even 
ones:  the  odd  clusters  combine  to  form  even  clusters. 

In  conclusion,  we  have  carried  out  self-consistent  ab  initio 
calculations  of  the  equilibrium  geometries,  energies  and  ionization 
potentials  of  small  Li  clusters  and  compared  the  results  with  the 
predictions  of  the  jellium  model.  We  find  that  clusters  with 
N  <  5  are  planar.  The  magnetism  of  these  clusters  is  not  governed 
by  the  exchange  interaction  alone.  The  N=4  and  N=5  clusters  in 
the  planar  structure  are  less  magnetic  than  in  the  three-dimen¬ 
sional  structure.  The  effect  of  correlation  has  been  confirmed 
as  crucial  in  determining  the  equilibrium  geometries  and  ioni¬ 
zation  potentials  of  small  metal  clusters.  Those  have  important 
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connection  with  the  interpretation  of  magic  numbers  observed  in 
the  experiments. 

This  work  was  supported  in  part  by  grants  from  Thomas  F. 
Jeffress  and  Kate  Miller  Jeffress  Trust,  National  Science  Foun¬ 
dation,  the  Petroleum  Research  Fund  administered  by  the  ACS,  and 
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Ficure  Caption 

Fig.  1.  Plot  of  second  derivative  of  E{N)  versus  N  1 
of  lithium  atoms.  E  is  given  in  atomic  Hart 

Fig.  2.  Shapes  of  optimized  geometries  of  lithium  at 

for  N  =  1-5. 
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